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We report values of R = a(e^ —* hadrons)/(T{e+e~ — > fi^/i^) for 85 center-of-mass energies 
between 2 and 5 GeV measured with the upgraded Beijing Spectrometer at the Beijing Electron- 
Positron Collider. 



In precision tests of the Standard Model (SM) Q], the 
quantities (^(M^), the QED running couphng constant 
evaluated at the Z pole, and = {g — 2)/2, the anoma- 
lous magnetic moment of the muon, are of fundamental 
importance. The dominant uncertainties in both a(M^) 
and a^^^ are due to the effects of hadronic vacuum po- 
larization, which cannot be reliably calculated in the low 
energy region. Instead, with the application of dispersion 
relations, experimentally measured R values are used to 
determine the vacuum polarization, where R is the lowest 
order cross section for e+e^ — s- 7* — > hadrons in units of 
the lowest-order QED cross section for e~^e~ ^J■^ : 
namely R = a{e^e^ — > hadrons )/cr(e^"e^ 



where a{e^ e M^M ) 



47ra2(0)/3s. 



Values of R in the center-of-mass (cm.) energy {Ecm) 
range below 5 GeV were measured about 20 years ago 
with a precision of 15 — 20% [^-|4|. In this paper, we 
report measurements of R at 85 cm. energies between 
2 and 4.8 GeV, with an average precision of 6.6% 
The measurements were carried out with the upgraded 
Beijing Spectrometer (BESII) at the Beijing Electron- 
Positron Collider (BEPC). 

Experimentally, the value of R is determined from the 
number of observed hadronic events, -/V/Jq^, by the rela- 



tion 



R 



luobs 



N, 



bg 



L ■ etrg ■ £had ■ {I + S) 



(1) 



where N^g is the number of beam-associated background 
events; Ei-^ii^ = r) are the numbers of lepton- 
pair events from one-photon processes and N^^ the num- 
ber of two-photon process events that are misidentified as 
hadronic events; L is the integrated luminosity; 6 is the 
effective initial state radiative (ISR) correction; ehad is 
the average detection efficiency for hadronic events; and 
etrg is the trigger efficiency. The triggers and the inte- 
grated luminosity measurement were the same as those 
used in a preliminary scan that measured i? at 6 energy 
points between 2.6 and 5 GeV Q. 

The hadronic event selection is similar with that used 
in the first R scan [Q but with improvements that in- 
clude: for good track selection, the distance of closest ap- 
proach requirement (< 18cm) of a track to the interaction 
point along the beam axis is not imposed; for event-level 
selection, the selected tracks must not all point into the 
forward {cos9 > 0) or the backward {cos 9 < 0) hemi- 
sphere. Some distributions comparing data and Monte 
Carlo data are shown in Figs. |^ (a)-(c). The cuts used for 
selecting hadronic events were varied over a wide range, 
e.g. |cos0| from 0.75 to 0.90, Esum from 0.24:Ebeam to 



0.32Ebeam, {Esum is the total dcpOSitcd energy, Etearn 

the beam energy) to estimate the systematic error arising 
from the event selection; this is the dominant component 
of the systematic error as indicated in Table |l|. 
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FIG. 1. Distributions for Ecm=3-0 GeV of (a) track mo- 
mentum; (b) track cos 6; (c) total energy deposited in the 
BSC; and (d) event vertex position along the beam (z) axis. 
Histograms and dots in (a)-(c) represent Monte Carlo and real 
data, respectively; the beam associated background in (c) has 
been removed by sideband subtraction. 

The numbers of hadronic events and beam-associated 
background events are determined by fitting the distri- 
bution of event vertices along the beam direction with 
a Gaussian to describe the hadronic events and a poly- 
nomial of degree one to three for the beam-associated 
background. This background varies from 3 to 10% of 
the selected hadronic event candidates, depending on the 
energy. The fit using a second degree polynomial, shown 
in Fig. |l| (d), turned out to be the best. The difference 
between using a polynomial of degree one or three to 
that of degree two is about 1%, which is included in the 
systematic error in the event selection. 

A special joint effort was made by the Lund group and 
the BES collaboration to develop the LUARLW gener- 
ator, which uses a formalism based on the Lund Model 
Area Law, but without the extreme-high-energy approx- 
imations used in JETSET's string fragmentation algo- 
rithm 1^. The final states simulated in LUARLW are 
exclusive, in contrast to JETSET, where they are inclu- 
sive. Above 3.77 GeV, the production of D, D* , Dg, and 
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D* is included in the generator according to the Eichten 
Model jo). A Monte Carlo event generator has been de- 
veloped to handle decays of the resonances in the radia- 
tive return processes e~^e~ "jJ/tp or ^^|J{2S) ]To[ |. 

The parameters in LUARLW are tuned to reproduce 
14 distributions of kinematic variables over the entire en- 
ergy region covered by the scan We find that one 
set of parameter values is required for the cm. energy 
region below open charm threshold, and that a second 
set is required for higher energies. In an alternative ap- 
proach, the parameter values were tuned point-by-point 
throughout the entire energy range. The detection ef- 
ficiencies determined using individually tuned parame- 
ters are consistent with those determined with globally 
tuned parameters to within 2%. This difference is in- 
cluded in the systematic errors. The detection efhcien- 
cies were also determined using JETSET74 for the ener- 
gies above 3 GeV. The difference between the JETSET74 
and LUARLW results is about 1%, and is also taken into 
account in estimating the systematic uncertainty. Fig- 
ure H (a) shows the variation of the detection efficiency 
as a function of cm. energy. 

We changed the fractions of D, D*, Dg, and D* pro- 
duction by 50% and find that the detection efficiency 
varies less than 1%. We also varied the fraction of the 
continuum under the broad resonances by 20%, and find 
the change of the detection efficiency is about 1%. These 
variations are included in the systematic errors. 
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FIG. 2. (a) The cm. energy dependence of the detection 
efficiency for hadronic events estimated using the LUARLW 
generator. The error bars are the total systematic errors, (b) 
The calculated radiative correction, and (c) the product of 
(a) and (b). 

Different schemes for the initial state radiative cor- 
rections were compared |T^-[T5[, as reported in ref. 0. 
Below charm threshold, the four different schemes agree 
with each other to within 1%, while above charm thresh- 
old, where resonances are important, the agreement is 



within 1 to 3%. The radiative correction used in this 
analysis is based on ref. and the differences with 

the other schemes are included in the systematic error 
]l6[ . In practice, the radiative effects in the detection 
efficiency were moved into radiative correction factor by 
making the replacement ehadi^ + S) ^ £(0)(1 + Sobs), 
where e{k) is the efficiency for events with a radiative 
photon of energy fc, and Sots contains a modification of 
the bremsstrahlung term to reflect the A:-dependencc of 
the hadronic acceptance. 

To calculate Sobsi a cutoff in s', the effective cm. en- 
ergy after ISR to produce hadrons, has to be made. 
In our calculation, the mimimum value of s' should be 
the threshold for producing two pions, corresponding to 
kmax = 1 - s'/s = (0.9805 - 0.9969) in the 2-5 GeV 
range. Our criteria to select hadronic events is such that 
e approaches zero when k is close to 0.90, which makes 
us insensitive to events with high ISR photon energy. 

In calculating the radiative correction for the narrow 
resonances J/ip and ^p{2S), the theoretical cross section 
is convoluted with the energy distribution of the collid- 
ing beams, which is treated as a Gaussian with a rel- 
ative beam energy spread of 1.32 x 10~^ Ecm {Ecm in 
GeV). For the broad resonances at 3770, 4040, 4160, and 
4416 MeV, the interferences and the energy-dependence 
of total widths were taken into consideration. Initially 
the resonance parameters from PDG2000 |l7j were used; 
then the parameters were allowed to vary and were de- 
termined from our fit. The calculation converged after a 
few iterations. 

We varied the input parameters (masses and widths) 
of the J/ip, ip{2S), and the broad resonances used in the 
radiative correction determination by one standard devi- 
ation from the values quoted in ref. , and find that the 
changes in the R value are less than 1% for most points. 
Points close to the resonance at 4.0 GeV have errors from 
1 to 1.7%. Figure H (b) shows the radiative correction as 
a function of cm. energy, where the structure at higher 
energy is related to the radiative tail of the ip(2S) and 
the broad resonances in this energy region. Tables ^ and 
^ list some of the values used in the determination of R 
and the contributions to the uncertainty in the value of 
i? at a few typical energy points in the scanned energy 
range, respectively. 

TABLE L Some values used in the determination of R at 
a few typical energy points. 





J^obs 
"had 


Na + 


L 


e(0) 


1 + Sobs 


R 


Stat. 


Sys. 


(GeV) 






(nb-i) 


(%) 






error 


error 


2.000 


1155.4 


19.5 


47.3 


49.50 


1.024 


2.18 


0.07 


0.18 


3.000 


2055.4 


24.3 


135.9 


67.55 


1.038 


2.21 


0.05 


0.11 


4.000 


768.7 


58.0 


48.9 


80.34 


1.055 


3.16 


0.14 


0.15 


4.800 


1215.3 


92.6 


84.4 


86.79 


1.113 


3.66 


0.14 


0.19 



Table III lists the values of R from this experiment. 
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TABLE II. Contributions to systematic errors: experimen- 
tal selection of hadronic events, luminosity determination, 
theoretical modeling of hadronic events, trigger efficiency, ra- 
diative corrections and total systematic error. All errors are 
in percentages (%). 





hadron 


L 


M.C. 


trigger 


radiative 


total 


(GeV) 


selection 




modeling 




correction 




2.000 


7.07 


2.81 


2.62 


0.5 


1.06 


8.13 


3.000 


3.30 


2.30 


2.66 


0.5 


1.32 


5.02 


4.000 


2.64 


2.43 


2.25 


0.5 


1.82 


4.64 


4.800 


3.58 


1.74 


3.05 


0.5 


1.02 


5.14 



They are displayed in Fig. |^, together with BESII val- 
ues from ref. 0| and those measured by MarkI, 772, and 
Pluto [|-|. The R values from BESII have an aver- 
age uncertainty of about 6.6%, which represents a factor 
of two to three improvement in precision in the 2 to 5 
GeV energy region. Of this error, 3.3% is common to 
all points. These improved measurements have a signifi- 
cant impact on the global fit to the electroweak data and 
the determination of the SM prediction for the mass of 
the Higgs particle . In addition, they are expected to 
provide an improvement in the precision of the calculated 
and test the QCD sum rules down 

to 2 GeV |21|,|2S 




> 
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FIG. 3. (a) A compilation of measurements of R in the 
cm. energy range from 1.4 to 5 GeV. (b) R values from this 
experiment in the resonance region between 3.7 and 4.6 GeV. 
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TABLE III. Values of R from this experiment; the first error is statistical, the second systematic {Ecm in GeV). 



E 

-'-'cm 


R 


E 


R 




R 




R 


2.000 


2.18 ±0.07 ±0.18 


3.890 


2.64 ±0.11 ±0.15 


4.120 


4.11 ± 0.24 ± 0.23 


4.340 


3.27 ±0.15 ±0.18 


2.200 


2.38 ±0.07 ±0.17 


3.930 


3.18 ±0.14 ±0.17 


4.130 


3.99 ±0.15 ±0.17 


4.350 


3.49 ±0.14 ±0.14 


2.400 


2.38 ±0.07 ±0.14 


3.940 


2.94 ±0.13 ±0.19 


4.140 


3.83 ±0.15 ±0.18 


4.360 


3.47 ±0.13 ±0.18 


2.500 


2.39 ±0.08 ±0.15 


3.950 


2.97 ±0.13 ±0.17 


4.150 


4.21 ±0.18 ±0.19 


4.380 


3.50 ±0.15 ±0.17 


2.600 


2.38 ±0.06 ±0.15 


3.960 


2.79 ±0.12 ±0.17 


4.160 


4.12 ±0.15 ±0.16 


4.390 


3.48 ±0.16 ±0.16 


2.700 


2.30 ±0.07 ±0.13 


3.970 


3.29 ±0.13 ±0.13 


4.170 


4.12 ±0.15 ±0.19 


4.400 


3.91 ±0.16 ±0.19 


2.800 


2.17 ±0.06 ±0.14 


3.980 


3.13 ±0.14 ±0.16 


4.180 


4.18 ±0.17 ±0.18 


4.410 


3.79 ±0.15 ±0.20 


2.900 


2.22 ±0.07 ±0.13 


3.990 


3.06 ±0.15 ±0.18 


4.190 


4.01 ±0.14 ±0.14 


4.420 


3.68 ±0.14 ±0.17 


3.000 


2.21 ±0.05 ±0.11 


4.000 


3.16 ±0.14 ±0.15 


4.200 


3.87 ±0.16 ±0.16 


4.430 


4.02 ±0.16 ±0.20 


3.700 


2.23 ± 0.08 ± 0.08 


4.010 


3.53 ±0.16 ±0.20 


4.210 


3.20 ±0.16 ±0.17 


4.440 


3.85 ±0.17 ±0.17 


3.730 


2.10 ±0.08 ±0.14 


4.020 


4.43 ± 0.16 ± 0.21 


4.220 


3.62 ±0.15 ±0.20 


4.450 


3.75 ±0.15 ±0.17 


3.750 


2.47 ±0.09 ±0.12 


4.027 


4.58 ±0.18 ±0.21 


4.230 


3.21 ±0.13 ±0.15 


4.460 


3.66 ±0.17 ±0.16 


3.760 


2.77 ±0.11 ±0.13 


4.030 


4.58 ± 0.20 ± 0.23 


4.240 


3.24 ±0.12 ±0.15 


4.480 


3.54 ±0.17 ±0.18 


3.764 


3.29 ±0.27 ±0.29 


4.033 


4.32 ±0.17 ±0.22 


4.245 


2.97 ±0.11 ±0.14 


4.500 


3.49 ±0.14 ±0.15 


3.768 


3.80 ± 0.33 ± 0.25 


4.040 


4.40 ±0.17 ±0.19 


4.250 


2.71 ±0.12 ±0.13 


4.520 


3.25 ±0.13 ±0.15 


3.770 


3.55 ±0.14 ±0.19 


4.050 


4.23 ±0.17 ±0.22 


4.255 


2.88 ±0.11 ±0.14 


4.540 


3.23 ±0.14 ±0.18 


3.772 


3.12 ±0.24 ±0.23 


4.060 


4.65 ±0.19 ±0.19 


4.260 


2.97 ±0.11 ±0.14 


4.560 


3.62 ±0.13 ±0.16 


3.776 


3.26 ±0.26 ±0.19 


4.070 


4.14 ±0.20 ±0.19 


4.265 


3.04 ± 0.13 ± 0.14 


4.600 


3.31 ±0.11 ±0.16 


3.780 


3.28 ±0.12 ±0.12 


4.080 


4.24 ±0.21 ±0.18 


4.270 


3.26 ± 0.12 ± 0.16 


4.800 


3.66 ±0.14 ±0.19 


3.790 


2.62 ±0.11 ±0.10 


4.090 


4.06 ±0.17 ±0.18 


4.280 


3.08 ±0.12 ±0.15 






3.810 


2.38 ±0.10 ±0.12 


4.100 


3.97 ±0.16 ±0.18 


4.300 


3.11 ±0.12 ±0.12 






3.850 


2.47 ±0.11 ±0.13 


4.110 


3.92 ±0.16 ±0.19 


4.320 


2.96 ±0.12 ±0.14 
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